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Abstract 

We study the decay width and CP asymmetry of the inclusive process b ^ sg {g de- 
notes gluon) in the multi Higgs doublet models with complex Yukawa couplings, including 
next to leading QCD corrections. We analyse the dependencies of the decay width and 
CP asymmetry on the scale /i and CP violating parameter 6. We observe that there exist 
an enhancement in the decay width and CP asymmetry is at the order of 10~^. 
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1 Introduction 



Rare B decays are induced by flavor changing neutral currents (FCNC) at loop level. Therefore 
they are phenomenologically rich and provide a comprehensive information about the theoretical 
models and the existing free parameters. The forthcoming experiments at SLAC, KEK B- 
factories, HERA-B and possible future accelerators stimulate the study of such decays since the 
large number of events can take place and various branching ratios, CP-violating asymmetries, 
polarization effects, etc., can be measured [|I|, 0]. 

Among B decay modes, inclusive b —>■ sg is interesting since it is theoretically clean and 
sensitive to new physics beyond the SM, like two Higgs doublet model (2HDM) P], minimal 
supersymmetric Standard model (MSSM) ^, etc. 

There are various studies on this process in the literature. The Branching ratio {Br) of 
b — > sg decay in the SM is Br{b sg) ~ 0.2% for on-shell gluon |^. This ratio can be 
enhanced with the addition of QCD corrections or by taking into account the extensions of 
the SM. The enhanced Br{b —>■ sg) is among the possible explanations for the semileptonic 
branching ratio Bsl and the average charm multiplicity . The theoretical predictions of Bsl 
are slightly different than the experimental measurements obtained at the T(4S') and Z° 
resonance Further the measured charm multiplicity rjc is smaller than the theoretical result. 
The enhancement of Br{B —>■ Xnocharm) and therefore Br{b sg) rate would explain the 
missing charm and Bsl problem [§. Further, Br[B —>■ rj'Xg) reported by CLEO JlOl stimulates 
to study on the enhancement of Br{b sg). 



In |Tl|, [T2[, the enhancement of Br {b sg) was obtained less than one order compared 



to the SM case in the framework of the 2HDM (Model I and II) for m/^± ~ 200 6*61^ and 



tan (3 ~ 5. The possibility of large Br in the supersymmetric models was studied in |T3[. In 
H Br was calculated in the model III and the prediction of the enhancement, at least one 
order larger compared to the SM one, makes it possible to describe the results coming from 
experiments [0]. In the case of time- like gluon, namely b sg* decay, Br should be consistent 



with the CLEO data [15 



Br {b sg*) < 6.8 % (1) 
and in IM , it was showed that the model III enhancement was not contradict with this data for 



light-like gluon case. The calculation of Br [b — » sg) with the addition of next to leading loga- 
rithmic (NLL) QCD corrections was done in and it was observed that this ratio enhanced 
by more than a factor of 2. 
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CP violating asymmetry (Acp) is another physical parameter which can give strong clues 
for the physics beyond the SM. The source of CP violating effects in the SM are complex 
Cabbibo-Cobayashi-Maskawa (CKM) matrix elements. Acp for the inclusive b ^ sg decay 
vanishes in the SM and this forces one to go beyond the SM to check if a measurable Acp is 
obtained. 

In this work, we study the decay width F and Acp oih ^ sg decay in the 3HDM and model 
III version of 2HDM. In these models, it is possible to enhance F and to get a measurable 
Acp- Since the Yukawa couplings for new physics can be chosen complex and the addition of 



NLL corrections brings additional complex quantities into the amplitude, theoretically, it 
is possible to get a considerable Acp, at the order of the magnitude 2 %. This effect is due to 
new physics beyond the SM, 3HDM and model III in our case. 

The paper is organized as follows: In Section 2, we give a brief summary of the model III 
and 3HDM{02) and present the expressions appearing in the calculation of the decay width 
of the inclusive b sg decay. Further we calculate the CP asymmetry Aqp of the process. 
Section 3 is devoted to discussion and our conclusions. 

2 The inclusive process 6 ^ in the framework of the 
multi Higgs doublet models 

In this section, we study NLL corrected b —* sg decay width and the CP violating effects in the 
framework of the multi Higgs doublet models (model III version of 2HDM and 3HDM) 

In the SM and model I and II 2HDM, the flavour changing neutral current at tree level is 
forbidden. However, they are permitted in the general 2HDM, so called model HI with new 
parameters, i.e. Yukawa couplings. The Yukawa interaction in this general case reads as 

where L and R denote chiral projections L{R) = 1/2(1 =1= 75), 0^, for k = 1,2, are the two 
scalar doublets, Qn are quark doublets, UjR and Djr are quark singlets, r]'lj'^ and ^fj^ are the 
matrices of the Yukawa couplings. The Flavor changing (FC) part of the interaction is given 
by 

^Y,FC = i^j^QiL4>2UjR + ifjQiL<p2DjR + h.C. . (3) 



The choice of 0i and 02 



\ 

v + Ho ^[ 



+ 



(4) 



and the vacuum expectation values, 

< 01 >= 



<(j)2>=0 



(5) 



allows us to carry the information about new physics in the doublet 02- Further, we take Hi, 
H2 as the mass eigenstates /iq, ^0 respectively. Note that, at tree level, there is no mixing 
among CP even neutral Higgs particles, namely the SM one, Hq, and beyond, /lo- 
in eq.(^) the couplings ^^'^ for the FC charged interactions are 



u 

ch 

D 

ch 



eutral 



VcKM ^neutral 



where ^^eSraZ is defined by the expression 



^U{D) _ (.rU{D).-i^u,{D)\/V(D) 
— \ ^R{L) ) S ^L(R) ■ 



L{R) 



(6) 



(7) 



where ^n^utrai is denoted as Here the charged couplings are the linear combinations of 

neutral couplings multiplied by Vckm matrix elements (see pO| for details). In the case of the 
general 3HDM, there is an additional Higgs doublet, ^3, and the Yukawa interaction can be 
written as 



+ pZ''QiL4>3UjR + pfiQiL(p3DjR + h.c. , 



(8) 



where pfj^ is the new Yukawa matrix having complex entries, in general. The similar choice of 
Higgs doublets 



1 

V2 







(9) 



with the vacuum expectation values. 



< 01 > = 







< 



_1_ ( 



>= ; < 0, >= 



(10) 



can be done and the information about new physics is carried beyond the SM in the last two 
doublets, 02 and 03. Further, we take Hi, H2, H3 and if4 as the mass eigenstates ho, Aq, h^, A'q 
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where /iq, A'q are new neutral Higgs bosons due to the additional Higgs doublet in the 3HDM 
(see 



The Yukawa interaction for the Flavor Changing (FC) part is 

^Y,FC = ^f,^Q^L<p2U,R + ^^,Q^L(p2D,R + pf,^ Q^L^P^U^n + pf,Q^L(t^zD ,R + h.C. , (11) 

where the charged couplings and p|^^ are 

Pch = Pn Vckm , 

Pch = ^CKM pN , (12) 



and 

C:U{D) _ (.rU{D).-i^u,{D)\/V{D) 
— \ ^R(L) ) ^ ^L{R) 1 

pT' = K[frY^^''%'^^ . (13) 

Since there exist additional charged Higgs particles, F^, and neutral Higgs bosons A'^ in 
the 3HDM, we introduce a new global 0(2) symmetry in the Higgs sector, considering three 
Higgs scalars as orthogonal vectors in a new space, which we call Higgs flavor space and we 
denote the Higgs flavor index by "m", where m = 1, 2, 3. The transformation reads 

02 = COS a 02 + sin a 03 , 



-sin a 02 + cos a 03 , (14) 



where a is the global parameter, which represents a rotation of the vectors 02 and 03 along 
the axis that 0i lies, in the Higgs flavor space. This symmetry ensures that the new particles 



are mass degenerate with their counterparts existing in model III (see for details). Further 
the Yukawa Lagrangian (eq.(§)) is invariant under this transformation if the Yukawa matrices 
satisfy the expressions 

F'U(D) cU{D) . -U{D) ■ 

iij = iij ' COS a + Pi/ 'sma, 

-iU{D) ^U{D) . . -U{D) /.p-x 

Pij = —qij stna + Pij cos a. [lb) 

and we get 



Therefore, it possible to parametrize the Yukawa matrices ^^(-^^ and p^(^) as : 

= lysine , 

p^ = ie'^sine, (17) 
where e'^(^) are real matrices satisfy the equation 

and the angle 9 is the source of CP violation. Here X^(^) = y^^X^^^) with X = ^, p, e 
and T denotes transpose operation. In eq. ([T7|) , we take complex to carry all CP violating 
effects in the third Higgs scalar. 

Now, we would like to continue the study of the inclusive process h sg. Our starting 



point is the recent calculation of NLL corrected decay width [|16 

T{h sg) = T^ + T^^""' , (19) 



where 



with 



Ci\D\\ (20) 



D = c's''ff + ^{cl'^ff-l^c',''ff + cUhln— + n) 
An 3 /i 



and r'"'*^'"'^ is the result for the finite part of bremsstrahlung corrections 



+ C',{kln^ + r2) + C's'"{{ls + 8 + Po)ln^ + rs)}, (21) 
p p 



where 



^brems ^ J ^j^^^j^^ ^ ^+ ^ ^-^ _^ ^+ _^ ^+ ^ ^+ ^ ^- ) ^ ^22) 



r+ = 48 Cl I Az23r '^fe K - 2 ^r) , 
r^t = yC2^|A^23pm^(m2-2E,E,) , 

= 24 (7^ I Aiirl^ rrib (16 mb - 16 Er - 8ml Eg + 6mb Eg Er + ml) , 



|2 2/2 
■^Oe//* A • 1 ^2 



r+ = 32 Cl C2 I Ai23 1"" ml {ml - 2 Eg Er) , 



r+ = 25Q Cl Re[Cr' '*At23] mi Eg Er , 
^C2Re[C!^^f*A^2s]mlEgEr 



28 ^ 

0,eff* 



T-^ = -96 C2 Re[C^8'"^ ^ * ^H7] mtimb (Eg + Er)~ 2 {E'^+E^ + Eg Er) 

+ 4^i^^i^^±^]/(E,E.) (23) 



Here = |C? and = Cl - \Cl, and Ci = ^\Gp V^.V^W^ and c<, = """I^^E'"' - (see m 



for details). In eqs. (pT]) and (^Sf ) the Wilson coefficients C^''^-^^ and (7^2) (eq- (p4|)) includes 
LL corrections and new physics effects enter into the expressions through the coefficients C^'^^^ 



2^1 "'^^ ~ ^2 6^1' "^^^ ~ 247r4 "tb-'tsl "^^^ "-^ " 96 647r2 

1(2) 

'8 

and Cl'^^^ (see eq. (pOD). The symbol is defined as 77 = "q'^^'' and (3q = 23/3. The 
vectors a^, /i^, e'j, /j', A;-, I'-, a'^, appearing during QCD corrections, and the Wilson coefficients 
Cl'^^i^mw), Ci'^^^mw) and Cg'^-^-^ {niw), the functions Ain and A«23 in eqs. (p3D, ri, r2, rg 
and the numbers li, I2, h in eq. (pTf) are given in . 



Now, we would like to start with the calculation of CP asymmetry for the inclusive decay 
underconsideration. The possible sources of CP violation in the model III (3HDM) are the 
complex Yukawa couplings. Our procedure is to neglect all Yukawa couplings except ^^^j and 
^N,bb (^TV.tt ^N,bb) (^ee eqs. ( p!7| , |T8|) and Discussion section) in the model III {3HDM{02)). 
Therefore, in the model III {3HDM{02)), only the combination i^^u^M,bb i^N,tt^N,bb) is respon- 
sible for Ac p. Using the definition of Acp 

A _ ^ sg) — T{b ^ sg) 

we get 

_|_ Qbr 

Acp = /m[e^,,,] j^:^ , (25) 

in the model III where and A^*^**") are the contributions coming from D-part (bremsstrahlung 

part) and they read as 

= ^ciAjImlA^] , 

n'"' = 2 C2 y dEgdEr{B5 /m[A23] + fie Im[Ah-j]) , 

A^ = 2c,{\A,\' + \^§^,,\'\A,\' + 2ArRe[^§^,,]Re[A,]} , 

A'^ = 2 C2 1 dE, dEr{B^ + Re[^§^J {B, ReiAt^s] + B, Re[At,j])} . (26) 

The functions ^5,6, 7 and -04,5,6 are defined as 

A, = {C'M [h + ln[^] + ri] + C'M [k + ln[^] + r^]) , 

/i fi 

A, = {v''/''A + A,) + ^[A, + xA,-%''/'''A, + A, 

47r 3 



+ (r/^^/23 ^ [(^^ + 8 + P0) H—] + rs] + A 



5J , 



A, = r^i^/^3A2{l + ^[r^^/^\-^ + (/8 + 8 + /?o)M^]+r8]}, (27) 

47r 3 /i 
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and 

54 = 5i + B2 (r/^^/'3 + A3) Re[A Z23] + 53 (7^1^/23^^ ^ /?e[A^i7] , 
Be = Bsv''^''A,, 



-Bi,2,3 appearing in eq. (pSf ) read 



52 = 32mlE,Er[8C, + ^C2] 



B 
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'28 



il7\ 



Here we use the parametrizations 



with 



1=1 



i=l 



and the Wilson coefficients 



4(x-l)4 8(x-l)= 



with 



mi 
mtrrib 



^^^^'^ = 2%. -1)3 ^A(^^'^'^ 

^ / N yt{.yt - 3) i/f 1 



The LL corrected Wilson coefficients C° and are 



C'M = {lv'^'' + lv-''^'')C'2iniw), (34) 



and 



C',{mw) = 1, 

C^{mw) = 0. (35) 
In eq. (pT]) the parameter x is given by 

X = 6.7441 (7^37/23 _ ^14/23) _ ^gg) 

In our calculations we take only CN,bb complex, ^^ ^^^ = e'^, where ^ is the CP violating 

parameter which is restricted by the experimental upper limit of the neutron electric dipole 
moment eq. (0). For 3HDM{02), it is necessary to make the following replacements: 









Im[^N,bb] - 


^N,bb ^'^'^ ^ 1 




Re[^N,bb] - 


-> e^^f,^ cos^ e , 




\fD |2 _ 
\s,N,bb\ 




(37) 



3 Discussion 

The general 3IIDM model contains large number of free parameters, such as masses of charged 
and neutral Higgs bosons, complex Yukawa matrices, i^j'^ , pif with quark family indices j. 
First, a new global 0(2) symmetry is introduced in the Higgs flavor space to connect the Yukawa 
matrices in the second and third doublet and to keep the masses of new charged (neutral) Higgs 
particles in the third doublet degenerate to the ones in the second doublet [|T^. Second, the 
Yukawa couplings, which are entries of Yukawa matrices, is restricted using the experimental 
measurements, namely, AF = 2 mixing, the p parameter [|TB| and the CLEO measurement [|T^], 



Br[B X,7) = (3.15 ± 0.35 ± 0.32) 10"^ . (38) 
The constraints for the FC couplings from AF = 2 processes and p parameter for the model 



III were investigated without QCD corrections [T^ and the following predictions are reached 



Kj = ^dj << 1 , id = 1,2,3 



where u{d) is up (down) quark and i,j are the generation numbers and further 

Xbb , Kb » 1 and Xtt , << 1 . (39) 
In the analysis, the ansatz proposed by Cheng and Sher, 

is used. Respecting these constraints and using the measurement by the CLEO [1^ Collabo- 
ration we neglect all Yukawa couplings except C,Ntt> ^Nbb ^^e model III. In 3HDM{02), the 
same restrictions are done by taking into account only the couplings e%^^ and e^f,b- 

This section is devoted to the study of the CP parameter sin6 and the scale fi dependencies 
of the decay width F and CP asymmetry of Aqp for the inclusive decay b sg, in the 
framework of the model III and 3HDM{02)- In our analysis, we restrict the parameters 6, 
^Yin and ^^^^ (CiVtt aiid e^^b) ^^e model III {3HDM{02)), using the constraint for {Cj'^-^l, 
0.257 < ICj^-^l < 0.439, coming from the CLEO data eq. (|3|) (see Q). Here C^^^ is the 
effective magnetic dipole type Wilson coefficient for b s'y vertex. The above restriction 
allows us to define a constraint region for the parameter (^tv,**) terms of ^§^i,b i^N,bb) 
and 6 in the the model III, {3HDM{02))- Further, in our numerical calculations we respect 
the constraint for the angle 6, due to the experimental upper limit of neutron electric dipole 
moment, namely 

dn < lO-^^e-cm (41) 

which places an upper bound on the couplings with the expression in the model III {3HDM (O2)): 
^MutN%) stnO < 1.0 (^(ef^,,6lv^,J stnH < 1.0) for m^^ ^ 200 GeV 0. 

Throughout these calculations, we take the charged Higgs mass mH± = 400 6*61^, and we 
use the input values given in Table ([^). 

Fig. |I] (^ is devoted to the sinO dependence of F for /i = mj,, i^^b = 40 {e^^b = 40 m^) 
and \rtb\ = \%^\ < 1 (|S^| < 1) in the model III {3HDM{02)). Here F is restricted between 

4jV,!)b '^N,bb 

solid (dashed) lines for Cj^^ > {Cj^^ < 0). As shown in Fig. |I|, the decay width F can 
reach (0.78 ± 0.06) x 10~^^ in the region 0.2 < sin e < 0.7 for Cf^ > and the possible 
enhancement, a factor of 4.2 compared to the SM one (0.185 ± 0.037) x 10~^^ GeV [0 can be 
reached. For 3HDM{02), the upper range for the decay width F is (0.79 ± 0.07) x 10^^*^ in 
the region 0.2 < sin 9 < 0.7 for Cj^-^ > and this leads to an enhancement, a factor of 4.3 
compared to the SM one. F decreases with increasing sin9 for Cj^^ > and it can get larger 
values compared to the C^^'^ < case, in both models. The sin 9 dependence of F is weak 
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Paramptpr 


Valiip 


rric 


1 /I ^'n^^^^ 




4.0 (^LreV ) 


It/ t/* I 


U.U4 




i ( ^^Vj^e V J 


mi^ 


80.26 (GeV) 


mz 


91.19 (GeV) 


Aqcd 


0.214 (GeV) 




0.117 



Table 1: The values of the input parameters used in the numerical calculations. 

for Cj^-^ < and for this case, it takes slightly smaller values in the 3HDM{02) compared 
to the ones in the model III. In our numerical calculations, we observe that the contribution 
of bremsstrahlung corrections are almost one order smaller as a magnitude compared to the 
rest. Further, the restriction regions for C^^^ > and C'^" < become more seperated with 
increasing values of the scale fi and this behaviour is strong in the 3HDM{02)- The scale 
dependence of F is weak for the values fi > 2 GeV and almost no dependence is observed for 
the large values of /i scale for both models, (see Figs. ^ and ^). 

In Fig. I and |], we present the sin 9 dependence of Acp for /i = m^, ^^bb = 40 mj, 
i^Nbb — 40 m;,) and \rtb\ < 1 in the model III {3HDM{02))- Here Acp is restricted in the 
region bounded by solid (dashed) lines for Cj^^ > {Cj^^ < 0). As shown in figures, \Acp\ 
reaches 2.5 % for sin 9 = 0.7 and all possible values of Acp are negative. However, for C^^^ < 0, 
the allowed region becomes broader and Acp can take both signs, even can vanish. For this 
case, |A(7p| reaches almost 1 % as an upper limit in both models. Further Acp is more sensitive 
to sin 9 in the 3HDM{02) compared to the model HI. 

Fig. 1^ and ^ represent the scale /i dependence of Acp for sin9 = 0.5, \^^bb\ i^Nbb) = 40 
and \rtb\ < 1 in both models underconsideration. The scale dependence of Acp is also weak for 
the values /i > 2 GeV similar to that of F. Here the increasing values of sin 9 cause to increase 
the size of restriction region. 

At this stage we give the numerical values of F and Acp for \^^bb\ = 40mfe (e^;,^ = 40 m^) 
and = m;, in the range 0.2 < sin9 < 0.7, for model HI [3HDM[02))'- 

0.72 (0.72) X 10"^^ GeV < F < 0.84 (0.86) x 10"^^ GeV (upper boundary) for Cj^^ > , 
0.28 (0.28) X 10"^^ GeV < F < 0.40 (0.42) x 10"^^ GeV (lower boundary) for Cj^^ > , 

F = 0.50 (0.48) X IQ-^^GeV (upper boundary) forGj^^ < , 
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r = 0.20 (0.20) X 10"^^ GeV (lower boundary) far Cj^^ < , 

(42) 



and 



0.0080 (0.0015) < \Acp\ < 0.0250 (0.0250) (upper boundary) for Cj^^ > , 
0.0050 (0.0010) < \Acp\ < 0.0170 (0.0165) (lower boundary) forC^^^ > , 
0.0020 (0.0010) < Acp < 0.0060 (0.0060) (upper boundary) far C^^^ < , 
-0.0100 (-0.0100) < Acp < -0.0020 (-0.0010) (lower boundary) farC^^^ < . 

(43) 

Now we would like to present our conclusions: 

• r can reach 0.84 (0.86) x 10~^^ in the model III {?,HDM{02)) and this is an enhancement 
a factor of 4 compared to the SM one. 

• A measurable CP asymmetry Aqp exists with the addition of NLL QCD corrections 
and choice of complex Yukawa couphng {pN,bh section 2)) in the model III 
{3HDM{02)). \Acp\ can be obtained at the order of the magnitude of %2.5. This 
physical parameter is coming from the new physics effects and it can give strong clues 
about the physics beyond the SM. 
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Figure 1: F as a function of sin9 for \rtb\ = ItS^I < 1, ^^^f, = 40mb and ^ = mi,. Here F is 

^N,bb ' 

restricted in the region bounded by solid (dashed) hues for C^-^^ > {C^^^ < 0), in the model 
in. Dotted hne represents the SM contribution. 



14 



0.9 



0.8 



1:::^ 0.7 

CD 



CO 

t 



o.c 



0.4 - 



X 

'o 0.3 
0.2 



T 1 1 1 1 1 1 1 r 



Q X I I I I I I I I I I 

■ 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 O.C 0.C5 0. 

sin9 

Figure 2: The same as Fig. | but for 3HDM{02). 
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Figure 3: The same as Fig. but F as a function of /i for sinO = 0.5. 
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Figure 7: The same as Fig. |^ but Acp as a function of fi. 
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